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ABSTRACT 

We discuss the predictability of a conservative system that drives a chaotic system 
with positive maximum Lyapunov exponent Ao, such as the erratic motion of an asteroid 
in the gravitational field of two bodies of much larger mass. We consider the case where 
in absence of feedback (restricted model), the driving system is regular and completely 
predictable. A small feedback of strength e still allows a good forecasting in the driving 
system up to a very long time T p ~ e _a , where a depends on the details of the system. 
The most interesting situation happens when the Lyapunov exponent of the total system 
is strongly chaotic with X tot ~ , practically independent of e. Therefore an exponential 
amplification of a small incertitude on the initial conditions in the driving system for any 
e ^ coexists with very long predictability times. The paradox stems from saturation 
effects in the evolution for the growth of the incertitude as illustrated in a simple model 
of coupled maps and in a system of three point vortices in a disk. 
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It is commonly believed that a sensible dependence on initial condition makes foreca- 
sting impossible even in systems with few degrees of freedom. This is the so-called butterfly 
effect discovered by Lorenz in a numerical simulation of a model of convection with three 
degrees of freedom [1]. Using his words, 'A butterfly moving its wings over Brazil might 
cause the formation of a tornado in Texas'. 

In general, a dynamical system is considered chaotic when there is an exponential 
amplification of an infinitesimal perturbation 80 on the initial conditions, with a mean 
time rate given by the inverse of the maximum Lyapunov exponent A [2]. Indeed, such a 
deterministic system is expected to be predictable on times t <C A -1 and to behave like 
a random system on larger times. The purpose of this letter is to show that there exists 
a wide class of dynamical systems where a large value of the Lyapunov exponent does 
not imply a short predictability time on a physically relevant part of the system. In this 
respect, one can speak of strong chaos A ^> without butterfly effect. 

In particular, we discuss the predictability of a conservative system that drives a stron- 
gly chaotic system with positive maximum Lyapunov exponent Ao- In absence of feedback 
the driving system is regular and completely predictable. A small feedback of strength e 
still allows to predict the future of the driving system up to a very long predictability time 
T p that diverges with e. The Lyapunov exponent of the total system is X tot ~ Ao, and so 
there is a regime of strong chaos for all e values. The absence of the butterfly effect stems 
from saturation effects in the evolution laws for the growth of an incertitude on the driven 
system 

To be explicit, let us consider a system with evolution given by two sets of equations 



with £ G R n and 17 G R m . The variables 17 are thus driven by a sub-system represented by 
the variables £, with a weak feedback of order e. A typical physical example is given by an 
asteroid moving in the gravitational field generated by two celestial bodies of much larger 
mass such as Jupiter and Sun [3]. Usually the feedback is neglected, and one considers 
the restricted three body problem, i.e. the 17 variables passively driven. However, a finer 
description should take into account even the influence of the asteroid on the evolution 
of the other two bodies, i.e. an 'active' driving where e / 0. This situation appears in 
many other phenomena, such as the active advection of a contaminant in a fluid, a simple 
example which will be discussed in this letter. 

The main properties of the system in absence of feedback, are the following: 

(1) the driver is an independent dynamical system that exhibits a regular evolution with 
zero Lyapunov exponent; 

(2) the driven sub-system is chaotic with positive maximum Lyapunov exponent, say Ao- 
In other terms, the behavior of the driver (£ variables) is completely predictable. 

However, as soon as e ^ 0, one should consider the total system which is obviously chaotic 
with a Lyapunov exponent X tot that, a part small correction of order e, is given by the 
Lyapunov exponent Aq of the chaotic driven sub-system. This means that there is an 





dr/ 
dt 



G(£, ri) 



(16) 



2 



exponential amplification of a small incertitude on the knowledge of the initial conditions 
even in the driver. This is an amazing result, since it is natural to expect that it is possible 
to forecast the behavior of the driver for very long times as e — > 0. Actually, intuition 
is correct while the Lyapunov analysis gives completely wrong hints on the predictability 
problem at difference with what commonly believed. The famous butterfly effect of Lorenz 
seems not to forbid the possibility of predicting the future of a part of the system. The 
paradox stems from saturation effects in the evolution for the growth of the incertitude. 
To fix notation and definitions, let us consider the evolution of the total system 

= /i(x) where x = (£, r,) G R n+m (2) 

The incertitude on its state is A(t) = x(t) — x'(t) where x and x' are trajectories starting 
from close initial conditions, i.e. |x(0) — x'(0)| = So- In the limit 80 — > 0, A can be confused 
with the the tangent vector z whose evolution equations are 



6?Z ^fi 

— = J(t) z where J t k(t) 



x(t) 



The maximum Lyapunov exponent is then defined as the exponential rate of the incertitude 
growth, 

A= lim lim - In I ) = lim - In |z(t)| (4) 

i->oo«5 ->0 t \ S J t^oo t 

It is worth stressing that the full equations for the evolution of an incertitude are non-linear: 

^ = J(t)A + 0(A 2 ) (5) 

so that the two limits in (4) cannot be interchanged. 

The predictability of the system is defined in terms of the allowed maximal ignorance 
on the state of the system, a tolerance parameter A max which must be fixed according the 
necessities of the observer. The predictability time is thus 

T p = sup {t such that|A(t')| < A max for t' < t} (6) 
t 

If A max <^ 1, (5) is well approximated by (3) and the predictability time can be roughly 
identified with the inverse Lyapunov exponent, since 

and the dependence on the initial error 80 and on the tolerance parameter A max is only 
logarithmic and can be safely ignored for many practical purposes. 

Suppose now to be interested only on the incertitude A^-* in the driver system. When 
e = 0, the Lyapunov exponent of the driver = so that it is fully predictable, (we 



3 



have in general Tp ~ S , the exponent (3 depending on the particular system), while 
the driven system has a Lyapunov exponent Ao > 0. However, for any e ^ the two 
sub-systems are coupled, and the global Lyapunov exponent is expected to be 

A = Ao + 0(e) (8) 

A direct application of (7) would give 

T«> ~ T p ~ J- (9) 

One thus obtains a singular limit, lim^o Tp (e) ^ Tp (e = 0). The troubles stem from 
the identification between incertitude A(t) and tangent vector z(t), which is not correct 
on long time scales. 

It is convenient to illustrate the problem in a simpler context as the main qualitative 
aspects of equations (la) and (lb) can be reproduced considering only the feedback effect 
in two coupled maps of the type 

6+i = LCt + eh{r ]t ) (10a) 

r/ f+ i = G{r] t ) (106) 

where the time t is an integer variable, £ £ R 2 , h = (^1,^2) is a vector, function of the 
variable rj £ R 1 whose evolution is ruled by a chaotic one-dimensional map G, and L is 
the linear operator corresponding to a rotation of an arbitrary angle 

_ / cos(8) -sin(8)\ 
~ \sin(6) cos(6) J ' 

When e = 0, one is left with two independent systems, one of them regular and of Hamil- 
tonian type, the other fully chaotic. 

These maps provide a simple, maybe the simplest, example of a system with two 
different temporal regimes: 

(A) short times where 8oexp(\o < 1 so that it is correct to ignore the nonlinear term in 
(5), so that A ~ z 

(B) long times where one should consider the full non-linear equation (5) for the incertitude 
growth. 

From the observer point of view, both these regimes might be interesting, according his 
particular exigence. If one is interested in forecasting the very fine details of the systems, 
the tolerance threshold A max could be quite small, hence T p ~ Aq~ . In general, however, a 
system is considered unpredictable when the incertitude is rather large (say discrimination 
between sun/rain in meteorology) and regime (B) is the relevant one. In that case, non- 
linear effects in (5) cannot be neglected and in order to give an analytic estimate of the 
predictability time we can use a stochastic model of the deterministic equations. Indeed 
the chaotic feedback on the evolution of the 'driver' system can be simulated by a random 
vector w, i.e. 

Zt +1 = LZ t + ew t (11) 
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The incertitude A\ is then given by the difference of two trajectories £ t and £' t originated 
by nearby initial conditions and evolves according the stochastic map 

A ( t % = LA^ + eW t (12) 

where Wt = w' t — w t . For short times t <C Xq 1 |ln$o| one cannot consider the 'noises' w t 
and w' t as uncorrelated so that the incertitude on the driver grows exponentially under 
the influence on the deterministic chaos given by the feedback, \A\ | ~ |A^| ~ exp(At), 
with A = Ao + 0(e). For long times t ^> Aq~ |ln$o|, the random variables are practically 
uncorrelated so that their difference Wt still acts as a noisy term. As a consequence the 
growth of the incertitude is diffusive, since the formal solution of (12) is 



A\ 6) = V + e J2 L ~ T W(r)j (13) 
and noting that L f is a unitary transformation, from (13) one can derive the bound 



et 1 ' 2 (14) 



where we have used the estimate | J^ =0 W T \ ~ t 1 ! 2 given by standard arguments borrowed 
from the central limit theorem. In conclusion, for our model maps (10a) and (10b), the 
predictability time on the driver diverges like 

T«> ~ e- 2 (15) 

although there is a regime of strong chaos since the total Lyapunov exponent A = Ao + 0(e) 
does not vanish with the strength of the feedback e. 

It is important to stress that the particular power of the diffusive law in a realistic 
model can be different from that of a random walk, since the deterministic chaos of the 
feedback could be better represented by random variables with appropriate correlations. 
The qualitative behavior exhibited by the stochastic model for the incertitude growth 
(exponential followed by a power law) can be tested in a direct numerical simulation of the 
coupled maps (10a) and (10b), where we choose the linear vector function for the feedback, 

M^) = {"Hi "n) (16a) 

and the logistic map at the Ulam point for the driving system, 

G(r]) = 4:r](l-r]) (166) 

with Lyapunov exponent Ao = In 2. Fig 1 shows the behaviors of the incertitude |A^| 
for the driver, at starting with an error So on the initial condition 770 and no error on the 
driver. At the beginning both |A^^| and |A^| grows exponentially. However, the phase 
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space available to the variable r/ is finite, so that |A^^| is bounded by a maximum value 
A M ~ 0(1). It will be attained at the time t = t* ~ A ln(AM/^o)- when the incertitude 
on the driver system is |A^| ~ eAm, and so much lower than the threshold. At larger 
times t > t* , the incertitude on the driven system remains practically constant and |A^| 
increases with a diffusive law of type (14) according to the mechanism described by the 
stochastic map (12). 

We have also studied a more realistic model of two coupled standard maps in action- 
angle variables / and 

/<;>=/<'> - £S in(C + ^) 

= 4" + 4!'. (I?) 

o(2) _ A2) , 7-(2) 

V t +i - Vt + h+i 

where K > e is a control parameter of order unity, such that the system (i^ 2 -*, #( 2 -*) is 
chaotic when e = 0. We do not discuss in details the results for these coupled maps, since 
they are qualitatively similar to those obtained for the simplified model (10). 

We now consider an application to a physical phenomenon, the motion of an ensemble 
of point vortices in a fluid. It is a classical problem in fluid mechanics, formally similar to 
the planetary motion in gravitational field. Both the systems are Hamiltonian with long 
range interactions. The main qualitative differences are that the Hamiltonian for point 
vortices does not contain a kinetic term and the motion is confined on the two-dimensional 
plane. The phase space for a collection of N vortices has thus 2N dimensions, related to 
the physical coordinates. 

Dynamical properties of point vortex systems have been studied by several authors 
interested in their chaotic motion and connection with two dimensional turbulence (see [4] 
for a review). The Hamiltonian theory for vortex motion inside a bounded domain was 
developed many years ago [5] for several boundaries. We are here interested in the motion 
in the unitary disk D for which the Hamiltonian takes the form 

N 

2^ 



47T 

i>j 



1 v-^ rf + rj — 2rirjCos6ij 1 ^ „ 

h = — Vr^-iog — — ^ — — + — Vrfiogu-r 



where r \ represent the circulation of the i-th vortex of coordinates x; = (x{ = ricosdi, m = 
ri-sinOi) and dij = 8i — 8j. The canonical conjugated variable are the scaled coordinates 
(TiXi, m) and the phase space is thus N times the configuration space D. The Hamiltonian 
(18) is invariant under rotations in the configuration space, thus the angular momentum 
is a second conserved quantity 

N 

L 2 = ^^ + yl) (is) 

i=l 

By general results of Hamiltonian mechanics, a system of two point vortices is always 
integrable, but we should expect chaotic motion for N > 2 vortices. 
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In the following we will consider N = 3 vortices, two of them carrying fixed circulation 
Ti = T2 = 1 and representing the driver, that without feedback is integrable. The third 
vortex, of circulation ^ = e represents the driven system which now makes the total system 
chaotic. In the limit e — > the third vortex becomes a passive particle (it is passively 
transported by the flow generated by the two unit vortices) and does not influence the 
integrable motion of the two vortices as for the three-body restricted problem in celestial 
mechanics. The restricted system is still chaotic, but the incertitude is confined to the 
passive tracer, while the motion of the two vortices is, in general, quasi-periodic. This 
limit is one of the simplest example of chaotic advection in two-dimensional flow and it 
will be studied in detail in another paper [6]. 

For our particular problem of three vortices the Hamiltonian can be rewritten in the 
following standard perturbation form 

H = # (xi,x 2 ) + e#i(xi,x 2 ,x 3 ) + e 2 # 2 (x 3 ) (20) 

The first term Ho describes the dynamics of the two unit vortices (and leads to integrable 
motion for e = 0); the term Hi represents the interaction with the small vortex of circu- 
lation e and the last term is due to the interaction of the third vortex with its own image. 
The 0(e) term is thus the perturbation to the integrable system Ho and we are reduced 
to the general framework described above if we identify £ = (xi,X2) and 17 = X3. The only 
difference is that now the 17 dynamics is not chaotic by itself, but chaoticity is induced by 
the interaction with the integrable system £. 

We now describe a typical simulation of error growth in the point vortex model which 
reproduces the effects obtained with the coupled maps model. We fix the value of the 
coupling constant (circulation of the third vortex) e = 10 -6 and the initial condition for 
the vortex positions are chosen in order to obtain chaotic motion with a global Lyapunov 
exponent A ~ 0.041. The initial incertitude on the coordinates of the small vortex is 
AM(0) = 10" 3 while we suppose to know the initial position of the two big vortices with 
a precision of A(«)(0) = 10" 8 . The saturat ion value for the incertitude is proportional to 
the disk radius, here Am ~ 1- 

We let the system evolve according to the Hamiltonian dynamics (18) for quite long 
time and we computed, at each time, the maximum value reached by the incertitude (we 
used the maximum because in this system incertitude show strong oscillations: this is 
a memory of the quasi-periodic behavior for e = 0). This represent the worst situation 
for making predictions. The upper scatter plot in fig 2 shows the time evolution of the 
incertitude for the driven system A^-'(t). We can recognize a short (t < 100) exponential 
growth until the nonlinear effect becomes important. At large time it > t* ~ 300) the 
incertitude saturates to its maximum value Am- The lower scatter plot represents the 
incertitude for the driver system of two vortices, A^\t). We can easily recognized the two 
expected limiting behavior represented by the two lines. For small times, the error grows 
exponentially, A^(t) ~ e ft where 7 ~ 0.064 is close to the global Lyapunov exponent. 
For long time, the power law behavior is recovered, A^\t) ~ et 1//a with a -1 ~ 0.88. 

In conclusion we must stress that all our results can generalized in a straightforward 
way to a weakly chaotic driver with a maximum Lyapunov exponent A^ <C Xq. In fact, 
the driver might be either conservative or dissipative. The important point is that the 
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dynamics of the driver has a much longer characteristic time than the driven system so that, 
for an observer interested in the predictability problem, the two systems can be practically 
decoupled. The Lyapunov analysis, although mathematically correct, does not capture the 
physically relevant features of the phenomenon, and the exponential dependence on initial 
conditions does not affect the possibility of forecasting the future of the driver on very long 
time scale. This is still true in systems with many different time scales instead of only two 
ones, as fully developed turbulence, where the inverse Lyapunov exponent is not related 
to the predictability time on the large length scale motion. 
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Figure captions 



Fig 1 Growth of the incertitude |A^| of the driver system in the coupled maps (10a) and 
(10b) as a function of time t , where the rotation angle 6 = 0.82099, the feedback 
strength e = 10 -5 , and the error on the initial condition of the driven system (10b) 
So = 10 -10 . Dashed line: exponential regime A^\t) = eSo exp(Ao t) where Ao = In 2. 
Full line: A^(t) = et 1 / 2 . 

Fig. 2 Incertitude growth for the point vortex model. Cross: maximum of the incertitude 
A^-'(t) on the third vortex. Diamond: maximum error A^-* on the driving system 
of two vortices. Dotted line: exponential regime A (i \t) ~ exp(7*) with 7 = 0.064. 
Dashed line: power law regime A^-* ~ et 1//a with a -1 = 0.88. 
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Figure 2 
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